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Abstract. Near normal incident far-infrared reflectivity spectra of [111] dysprosium 
titanate (Dy2Ti207) single crystal have been measured at different temperatures. 
Seven phonon modes (eight at low temperature) are identified at frequency below 1000 
cm~^. Optical conductivity spectra are obtained by fitting all the reflectivity spectra 
with the factorized form of the dielectric function. Both the Born effective charges 
and the static optical primitivity are found to increase with decreasing temperature. 
Moreover, phonon linewidth narrowering and phonon modes shift with decreasing 
temperature are also observed, which may result from enhanced charge localization. 
The redshift of several low frequency modes is attributed to the spin-phonon coupling. 
All observed optical properties can be explained within the framework of nearest 
neighbor ferromagnetic(FM) spin ice model. 
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1. Introduction 

Recently, there has been a surge of interest in the properties of pyrochlore compound 
Dy2Ti207, which is considered to be a model system of "spin ice" materials. "Spin 
ice" materials governed by the same statistical mechanics of so-called "ice rule" 
as the hydrogen atoms in the ground state of ordinary hexagonal ice have 
macroscopically degenerate ground states down to almost zero temperature. [H 12 El 
m El El 13 IHl Ej Experimentally, the observed value of (l/2)Rln(3/2) through specific 
heat measurement ^ is consistent with what is expected by pauling's theory, while the 
spin entropy only freezes out below about 4 K.|2j From magnetic susceptibility studies, 
a strongly frequency dependent cooperative spin freezing is observed at about 16 K, 
which is associated with a very narrow distribution of spin relaxation times and a sharp 
drop at about 2 K.^Ej Neutron scattering studies performed by Fennell et al.jlj also 
well demonstrate the spin ice state and the coexistence of long range ferromagnetic and 
short range antiferromagnetic order in a magnetic field applied along the [110] axis of 
Dy2Ti207. In theoretical aspect, R.G. Melko and coworkers [7j report numerical results 
on the low temperature properties of the dipolar spin ice model by the multicanonical 
Monte Carlo (MC) method and they find a first order transition to a long-range ordered 
phase. Other researchers [SJ also confirm the existence of the transition under a magnetic 
field along the [110] axis with MC simulation. 

Dy2Ti207 has a typical A2B2O7 structure with the space group (Fd3m, Oh), 
No. 227. B cation is sixfold coordinated and locates at the center of the distorted 
octahedron formed by corner O ions. The A-site ions, i.e., Dy'^+(a magnetic rare-earth 
ion with effective spin S=l/2) resides on a lattice of corner-sharing tetrahedral, as shown 
in figure 1 of reference 3. In the Dy'^+ sublattice topology the spin configuration with 
two spins pointing directly towards while two spins pointing directly away from the 
center of the tetrahedral corresponds to the analogous proton disorder in real ice. Then 
a ferromagnetic and dipolar nearest-neighbor spin-spin interaction leads to a strong 
geometrical frustration, preventing the system from long range ordering. As a result, 
high degeneracy of spin states will occur at low temperature. 

Owing to the Ising anisotropy arising from crystal-field effect in the pyrochlore 
lattice, the macroscopically degenerate states when the field is along the [111] direction 
are different from those for other field directions. For the [111] direction, the frustration 
structure changes from that of a three-dimensional pyrochlore to that of a two- 
dimensional Kagome-like lattice with constraint, leading to different values of the zero- 
point entropy. [5j The combination of ferromagnetic coupling and Ising anisotropy may 
be a reason for the strong frustration. But how the spins interact with each other is not 
well understood in the spin ice materials. 

Up to now, there has been no experimental report on the infrared optical properties 
of Dy2Ti207. IR spectroscopy can give insight to the dynamical processes related to 
phonon, charge carrier, and spin. Motivated by this situation, we studied the infrared 
reflectivity of Dy2Ti207 single crystal on the [111] plane from room temperature down 
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to 10 K. The temperature dependence of phonon modes is obtained, and the role played 
by spin-phonon coupling on the phonon modes is being discussed. 

2. Experimental details 

Dy2Ti207 single crystal was prepared by the floating zone method, using an infrared 
furnace equipped with two elliptical mirrors. Before the single crystal growth, the 
polycrystalline rods was prepared by a standard solid state reaction. Stoichiometric 
mixture of DysOg (CERAC, 99.99%) and TiOa (CERAC, 99.9%) was heated in air at 
1250°C for five days with intermediate regrinding to ensure complete reaction. The 
typical growth condition was 4.5 mm/h for both feed and growth speeds. To avoid 
oxygen deficiency, the single crystal was grown in an O2 atmosphere of 0.3 MPa. The 
single crystal obtained was translucent yellow. Powder x-ray diffraction measurements 
on the crystal confirmed that the product was a single phase with cubic pyrochlore 
structure and the parameter of the unit cell was a= IO.II22A. The principal axes were 
determined using Laue diffraction pattern. 

The crystal was cleaved and one-side polished with the surface parallel to the [111] 
plane. The reflectivity spectra R(a') at different temperatures were measured at near- 
normal incidence of about 8° on a Bomen DAS Fourier transform infrared spectrometer, 
in the range from 50 to 6500 cm~^. In the far infrared region, a He cooled bolometer and 
6 mm mylar beamsplitter were used. In the middle infrared region, a liquid nitrogen 
cooled mercury cadmium telluride detector and KBr beamsplitter were utilized. To 
obtain the absolute reflectivity, an evaporated golden mirror was served as a reference. 
Spectra were collected with a resolution of 4 cm~^. The samples were mounted in a 
continuous helium flow cryostat in which the temperature could be varied between 300 
and 10 K. 

3. Results and discussion 

The temperature dependent reflectivity of the Dy2Ti207 single crystal shown in figure 
1 is typical that of a nonmetallic system. The sharp features in the reflectivity spectra 
are due to the unscreened infrared active optical phonon modes, above the highest 
observed lattice vibration frequency, the reflectivity is flat and featureless up to the 
highest measured frequency. It can be seen that the reflectivity within the reststrahlen 
bands increases with decreasing temperature. At high temperature seven modes can 
easily be identified. 

The mode frequencies are obtained from a damped harmonic oscillator fit of 
reflectivity spectra with a complex dielectric function e{uj) in the factorized form 
(generalized Lyddane-Sachs- Teller relation) [TIH ITTj: 




up' — i'JjTO^ 



(1) 



Far-infrared optical properties of the pyrochlore spin ice compound Dy2Ti20j 



4 



where £00 represents the high frequency dielectric constant (at frequencies large 
compared with the lattice vibration frequencies but small compared with the 
electronic transition frequencies), VtjLo and VtjTo are the longitudinal and transverse 
eigenfrequencies of the jth optical phonon mode, 7^x0 and 7jto represent their 
longitudinal and transverse damping constants. Using this dielectric function, all the 
reflectivity spectra in our measurements were fitted with the well known Fresnel formula 
for reflectivity of a half space medium in vacuum [T^: 



Based on equation (1) and (2), a fit of R(w) to the observed reflectivity spectra can be 
obtained with a proper choice of the model parameters fijLO; ^jTO; IjLOi IjTO and 
. Adjustment of the parameters is made by trial and error fitting of formula (2) to the 
experimental spectra. This method yields not only e{u) but also the model parameters 
which characterize the infrared active phonons. 

From a group theoretical analvsis.|13j a compound with the pyrochlore structure 
symmetry should have the vibrational phonon modes at the Brillouin zone center F 
point with 



Among these 26 normal modes, only Aig, Eg, 4 are Raman active, 7Fiu are infrared 
active, and one Fi^ acoustical. Consequently, the observed phonons in the spectra of 
Dy2Ti207 at low temperature except that at about 610 cm~^ were denoted sequentially 
from low to high frequencies by Fi^^ to Fi„''. The number is in good agreement with 
the group theoretical analysis. The mode at 610 cm"^, denoted by Fi^'''*, is weak in 
intensity and does not change with the temperature. It is probably originated from two 
phonon absorption because of its low intensity. jH] 

The fitting results of the oscillator parameters for all the phonon modes of Dy2Ti207 
are summarized in Table ^ The high frequency dielectric constant is adopted as 5.0. 
Figure 2 shows the experimental and fitting curves at two representative temperatures of 
300 K and 10 K. The fitting results with the Lorentz oscillators could explain the global 
features of the phonon spectra reasonably well, even though some difference between 
the experimental and the fitted data remains. A small deviation of the calculated curves 
from the experimental one is seen on the low frequency edge of Fi^j^ mode in the high 
temperature. A possible reason is the reflection from the rear surface of the sample due 
to its transparency at low frequencies as well as the cutoff of the working frequency of 
the detector. However, the deviation has few impact on other phonon fitting parameters 
in other frequency range. 

The real part of the optical conductivity ai{uj) can be extracted from £2('^); where 
cji (a;) =1.6638 • 10~^a'£2('^), here u is in unit of cm~^ and cti in unit of cm~^, 
the corresponding result is shown in figure 3.^3] The conductivity is dominated by 
seven peak structures due to optical phonon absorption, no free carrier contribution 
can be observed. For normal materials, with decreasing T, the anharmonic thermal 




(2) 



r=8 Fi„+4 F2„+2 Fi, +4 F^g +3 E„ + + Ai,, +3 A2„. 
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motion would decrease, which results in a decrease in the lattice constant. Then 
the phonons should shift to higher frequencies and their linewidth should become 
narrower. The Fi„^ mode has clearly demonstrated this effect. Its center frequency 
shifts from 439 cm~^ to 452 cm~^ as the temperature decreases from 300 K down 
to 10 K. Moreover, the phonon peak is more and more prominent and separated 
from the nearby mode Fi„^ . But it is puzzling that other modes do not exhibit 
the same behavior, the three most prominent phonons, Fi„^, Yi^^ and Fi„^, show 
discernible redshifts. Coupling of phonons to magnetic excitations in solids may result 
in a change of phonon self-energy, i.e., the frequency due to spin-phonon interaction. 
Therefore, spin-phonon couplingjinj closely correlated with the Dy^"^ (its effective spin 
S=l/2) could be a possible origin for the abnormal temperature dependence of the 
phonon frequency. As we mentioned in the introduction, Dy^"*" cations can have 
appreciable magnetic moments. Consequently, ferromagnetic nearest neighbor dipole- 
dipole interactions will present in spin ice Dy2Ti207 in which large dipole interactions 
have been suggested to be responsible for the spin ice behavior. |17j These effects have 
been quantitatively confirmed by both experimental and theoretical approaches. [H] At 
the nearest neighboring sites, the exchange interactions between the magnetic atoms 
become the strongest. Subsequently, the crystal potential U can be described as 
following: 



where x is the atomic displacement from the equilibrium position in the oscillator 
model. In the second term, Jj^ is the exchange energy constant which is a function 
of the structural parameters, such as the bond lengths between Dy^+ cations and the 
corresponding bond angles mediated by the O ions. Sj represents the spin of Dy^"*" cation 
at the ith site. < SiSj > denotes a statistical average over the adjacent spins. The 
harmonic force constant derived from its second derivative formula reads: 



Note that the second term represents the spin-phonon coupling, which suggests that 
the phonon frequency should have an additional contribution. As for the spin-phonon 
coupling coefficient, Y.ij{~g^) can be different for each phonon and can have either a 
positive or a negative sign. Furthermore, different phonon frequencies will have redshift 
or blueshift in the optical conductivity spectra. As the temperature decreases, spin 
fluctuation becomes weaker and spin-phonon coupling stronger. It is important to 
notice that although ferromagnetic correlations are short ranged without long range 
order in spin ice compounds, from neutron scattering measurements Harris et al.[TH] 
have shown that the range of ferromagnetic order increases at lower temperature. As a 
result, spin-phonon coupling range is wider so that the phonon frequency shows redshift 
or blueshift, which is in agreement with what is observed in our optical conductivity 
spectra. 




(3) 




(4) 



Far-infrared optical properties of the pyrochlore spin ice compound Dy2Ti20j 



6 



Now we focus on the remarkable change of the spectra weight which is proportional 
to the area under the optical conductivity peak. It is well known that, by decreasing 
the temperature, thermal fluctuation reduces, which may result in narrowing in the 
linewidth and enhancement in the oscillator strength without any changes in the bonding 
or coordination. Almost all phonon modes exhibit the theoretically predicted behavior 
in our optical conductivity spectra. But the spectral weight should not change. The 
anomalous increase in oscillator strength of the low-frequency modes must be taken 
sufficient care of. Optical sum rules provide a powerful tool with which to analyze the 
behavior of free carriers and bound excitations. |19j The partial conductivity sum rules 
corresponding primarily to a single class of absorption such as excitation of phonons, 
conduction, valence, or core electrons have been developed. For oscillator states, the 
partial conductivity sum rule can be expressed as [20j 



effective plasma frequency associated with the isolated phonon absorption, respectively. 
In the above equation, a is in unit of Vf^cmT^. The integral region of the jth oscillator 
should be chosen to cover the full spectral weight. The dramatic increase in the spectra 
weight has implications for the distribution of charge and the change in local strength 
of the binding charge. The effective charge of Dy, Ti and O ions in the unit cell with k 
atoms can be determined through the following Equation (6) (in Gauss unit system) [21] 



where Vc is the unit cell volume, j and k index the lattice modes and the atoms with 
mass Mfc, respectively. For the effective charge, there is a general expression J2k^k~^- 
In Dy2Ti207, oxygen is the lightest element, therefore, in the right side of equation 
(6), the summation is dominated by O ion item and the terms for Dy ion and Ti ion may 
be neglected. The change in the effective charge is associated mainly with the oxygen 
(i.e., Z^RiZ*). Combining equation (5) and equation (6), we obtain the value for Z* as 
shown in figure 4. The absolute value of Z* is of less significance, what is important is the 
temperature dependence of the deduced value for Z*, which increases with decreasing 
temperature. When the incident hght couples to the induced dipole moments created by 
the atomic displacements associated with a normal mode, if the Born effective charge per 
oxygen atom Z* is increasing, then the size of the induced dipole moment and the optical 
absorption will also increase, which well explain the above temperature dependence of 
the optical conductivity. On the other hand, the increase in Z* means a change in the 
bond length and bond angle between O ion and cations so that the bond lengths between 
Dy^+ cations and the corresponding bond angles mediated by the O ions will change. 
Subsequently, Jjj and < SiSj > have different changes associated with phonon central 
frequency shift as we observed above. In summary, the change in the effective charge 
arising from electrical charge localization strengthening is the fundamental reason why 




(5) 



where Ua, ui^ and ujpj (in unit of cm ^) are the integral lower and upper limits and 




47r {Zlef 



(6) 
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the spectra weight increases and phonon frequency shifts. Prom our experimental results, 
within the framework of the ferromagnetic spin ice model, wc consider that the increase 
in Z* is closely related with nearest-neighbor ferromagnetic interaction strengthening 
and/or widening in the correlation range. In this simple physical picture, the increasing 
FM exchange between Dy^+ cations will lead to a decrease of the distance between 
mediated O ion and Dy^+, which enhance the increase in Z*, in agreement with our 
experimental observation. 

LST relationship also provides an effective examination on static optical 
permittivities £0 in the approximation of zero phonon frequency. When cu—O , the 
equation evolutes into the following form: 

too i ^HTO 

Prom Equation (7) Eq of Dy2Ti207 single crystal at different temperatures are also 
obtained, as shown in the inset of figure 4. It exhibits the same temperature dependence 
as Z*, which suggests that the temperature dependence of £0 can be originated from that 
of Z* and essentially, charge localization arising from PM exchange plays an important 
role in the temperature variation of both Eq and Z*. 

4. Conclusions 

In conclusion, PIR response of Dy2Ti207 single crystals has been studied under different 
temperatures. Seven infrared active phonons have been observed, which agrees well 
with a group analysis. All the spectra arc fitted with the oscillator model and an 
excellent agreement between experimental and calculated spectra is obtained. Of all 
phonons, one shows discernible blueshift and three obvious redshifts with decreasing 
temperature. The shifts are attributed to the spin-phonon coupling in a geometrically 
frustrated configuration in the spin ice material. The oscillation strength of the low 
frequency modes increase dramatically at low temperature, indicating that the Born 
effective charges are increasing in the unit cell. We propose that the similar temperature 
dependence of the static optical permittivity with Born effective charge may originate 
from the intrinsic charge localization resulted from the nearest-neighbor PM interaction. 
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Figure captions 

Figure 1 The temperature dependence of the reflectance of Dy2Ti207 single crystal 
in the range of 50-1000 cm^^. 

Figure 2 Representative experimental and fitted refiectivity spectra of Dy2Ti207 at 
(a) 300 K and (b) 10 K. Solid lines represent experimental data and dot lines are the 
fitting results. 

Figure 3 The real part of the temperature dependent optical conductivity of 
Dy2Ti207 in the far-infrared region. 

Figure 4 The temperature dependence of Born effective charge per oxygen atom in 
Dy2Ti207. The dots represent the deduced values at various temperatures between 300 
and 10 K. Inset: the temperature dependence of the static optical permittivities. 
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Table 1. The phonon parameters for the Lorentzian fits to the conductivity of 
Dy2Ti207 single crystal at different temperatures. All units are in cm~^. 
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Fig.l C.Z. Bietal. 




Fig.2(a) C.Z. Bi et al. 
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Fig.2(b) C.Z. Bi et al. 
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Fig.3 C.Z. Bi et al. 




Fig.4 C.Z. Bi et al. 



